BOUNDARY VALUE PROBLEMS ON PLANAR GRAPHS AND FLAT 
SURFACES WITH INTEGER CONE SINGULARITIES, I: THE 

DIRICHLET PROBLEM 

SA'AR HERSONSKY 

Abstract. Consider a planar, bounded, m-connected region fi, and let d£l be its boundary. 
Let T be a cellular decomposition of U <9fi, where each 2-cell is either a triangle or a 
quadrilateral. From these data and a conductance function we construct a canonical pair 
(5, /) where S is a genus (m — 1) singular flat surface tiled by rectangles and / is an energy 
preserving mapping from onto S. By a singular flat surface, we will mean a surface 
which carries a metric structure locally modeled on the Euclidean plane, except at a finite 
number of points. These points have cone singularities, and the cone angle is allowed to 
take any positive value (see for instance [28] for an excellent survey). Our realization may 
be considered as a discrete uniformization of planar bounded regions. 



0. Introduction 

Consider a planar, bounded, m-connected region fl, and let dfl be its boundary. Let T be 
a cellular decomposition of flUdfl, where each cell is either a triangle or a quadrilateral. Let 
dfl = Ex U E2, where E\ is the outermost component of dfl. Invoke a conductance function 
on making it a finite network, and use it to define a combintorial Laplacian A on T^°\ 
Let k be a positive constant. Let g be the solution of a Dirichlet boundary value problem 
(D-BVP) defined on T^' and determined by requiring that g\ El = k,g\ E2 = and Ag = 
at every interior vertex of T^°\ i.e. g is combinatorially harmonic. Furthermore, let E(g) 
be the Dirichlet energy of g. Following the notation in [T3], let ^{x) denote the normal 
derivative of g at the vertex x G dfl. Our first result is: 

Theorem 0.1. (Discrete uniformization of a pair of pants) Let V be a bounded, triply con- 
nected planar region and let (fl, dfl, T) = (V, dV = E x U E 2 , T) where E 2 = E\ U E\ . Let 
S-p be a singular flat pair of pants with exactly one singular point Q having An as its cone 
angle such that 

(1) Length Euclidean (S P ) Bl = J2 xe E 1 

(2) Length Euclidean (S P ) s i = - ExefiJ ff( x )> 

(3) Length Euclidean (,5'-p) S | = -E^ef f;( x )> and 

(4) the Euclidean length of a shortest geodesic connecting E\ to either E\ or E\ is k, 
where (S'p)^, (Sp) E i and (S-p) E 2 are the boundary components of S-p. Then there exists a 
mapping f which associates to each edge in a unique Euclidean rectangle in S-p, in such 
a way that the collection of these rectangles forms a tiling of S-p. Furthermore, f is energy 
preserving in the sense that E(g) = Area(Sp), and f is boundary preserving. 
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Throughout this paper a Euclidean rectangle will denote the image under an isometry of a 
planar Euclidean rectangle. For instance, most of the image rectangles that we will construct 
embed in a flat Euclidean cylinder. These will further be glued in a way that will not distort 
the Euclidean structure (see gand §@for the details). In our setting, boundary preserving 
means that the rectangle associated to an edge [u, v] with u G dQ has one of its edges on a 
corresponding boundary component of the singular surface. 

A singular flat, genus zero compact surface with m > 3 boundary components with conical 
singularities, will be called a ladder of singular pairs of pants. We now state the main result 
of this paper: 

Theorem 0.2. (The general case) Let (Q,dQ,T) be a bounded, m-connected, planar region 
with E-2 — E\ U E\ . . . U E™~ x . Let Sq be a ladder of singular pairs of pants such that 

(1) Length Euclidean (Sn) Bl = E*eEi ^( x )> and 

(2) Length EucMean (# n ) B i = - ^2 xeE i f£(a:), for i = 1, • • • , m - 1, 

where (Sq)e 1 and (Sn)^i, for i = 1, • • • , m — 1, are the boundary components of Sq. Then, 
there exists a mapping f which associates to each edge in a unique Euclidean rectangle 
in Sq i n such a way that the collection of these rectangles forms a tiling of Sq. Furthermore, 
f is boundary preserving, and f is energy preserving in the sense that E(g) = Area(Sn). 

The following Corollary is straightforward, thus establishing the statement in the abstract 
of this paper 



Corollary 0.3. Under the assumptions of Theorem 0.2, there exists a canonical pair (S,f), 
where S is a flat surface with conical singularities of genus (m — 1) tiled by Euclidean 
rectangles and f is an energy preserving mapping from T^ 1 ' into S , in the sense that 
2E(g) = Area(S'). Moreover, S admits a pair of pants decomposition whose dividing curves 



have Euclidean lengths given by (1) — (3) of Theorem 0.2 



Proof. Given (Q,dQ,T), glue together two copies of Sn (their existence is guaranteed by 



Theorem 0.2) along corresponding boundary components. This results in a flat surface 

S = Sq |J Sq of genus (m — 1) and a mapping / which restricts to / on each copy. 
an 

□ 



In the course of the proofs of Theorem 0.1 and Theorem 0.2, it will become apparent 



that the number of singular points and their cone angles, as well as the lengths of shortest 
geodesies between boundary curves in the ladder, may be explicitly determined. In partic- 
ular, the cone angles obtained by our construction are always even integer multiples of 7r 



(see Equation (4.12) for the actual computation). Some of these surfaces, those with trivial 
monodromy, are called translation surfaces and for excellent accounts see for instance }21j . 
[25] and [30]. Also, the dimensions of each rectangle are determined by the given D-BVP 
problem on T^. Concretely, for [u,v] G T^ l \ the associated rectangle will have its height 
equal to (g(u) — g(v)) and its width equal to c(u,v)(g(u) —g(v)), when g(u) > giy). Some of 
the rectangles constructed are not embedded. We will comment on this point (which is also 



transparent in the proofs) in Remark (4.15). In a snapshot, some of the rectangles which 
arise from intersection of edges with singular level curves are not embedded. 

The following theorem is foundational and serves as a building block in the proofs of all 
of the above theorems. We prove: 
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Theorem 0.4. (Discrete uniformization of an annulus [10J) Let A be an annulus and let 
(Q,dQ,T) = (A,dA = E\ U E2,T). Let be a straight Euclidean cylinder with height 
H = k and circumference 

Then there exists a mapping f which associates to each edge in a unique embedded 
Euclidean rectangle in Sjs, in such a way that the collection of these rectangles forms a tiling 
of Sa. Furthermore, f is boundary preserving, and f is energy preserving in the sense that 
E{g) = kre&{S A ). 

Given (Q,dQ,T), we will work with the natural afhne structure induced by the cellu- 
lar decomposition. Let us denote this cell complex endowed with this affine structure by 
C(Q,dQ,T). There is a common thread in our proofs of the theorems above. Extend g 
to an affine map g defined on T. This results in a piecewise linear structure on T. Next, 
study the level curves of g on a 2-dimensional complex which is homotopically equivalent 
to C(f2, <9f2, T), embedded in R 3 , obtained by using g as a height function on C(Cl,dCl,T). 
We will work with the level curves of g or equivalently, with their projection on C(Q, dQ, T). 
The topological structure of the level curves associated with the solutions will be carefully 



studied. For example, the level curves in the case of an annulus (Theorem 0.4) are simple, 
piecewise linear closed curves, all of which are in the (free) homotopy class determined by 
Ei. One nice consequence of this is that all the rectangles constructed in the proof of The- 
orem CL4 are embedded. In the proof of Theorem 0T, we will show that all the level curves 
associated to values in [0, ki], for some constant ki < k, are simple, piecewise linear closed 
curves in either the (free) homotopy class of E\ or the class determined by E\. However, 
for the value ki, it will be proved that the (unique) associated level curve is a figure eight. 
Furthermore, any level curve of g for a value which is larger than ki and smaller than or 
equal to A;, is a simple closed curve in the (free) homotopy class determined E\. The basic 



idea of Theorem 0.1 is to cut (Q,dQ,T) along the (projection of) a figure eight curve, tile 
each cylinder according to Theorem |0.4[ and glue back. 

We will often work with a series of modified boundary value problems. Each is a slight 
modification of the initial problem. This important feature of the proofs is due to the fact 
that the level sets of the original boundary value problem (defined on T^) intersect in 
a set which is much larger than (see §2 for details). 



In fact, once Theorem CL4 is proved, we proceed to prove Theorem (L2 by an inductive 
process. Unlike common proofs in the theory of surfaces, in which a surface is cut along closed 
geodesies yielding a pair of pants decomposition, we keep cutting our region along particular 
(projection of) singular level curves until we encounter a planar pair of pants or an annulus. 
This is a subtle point, arising from the fact that our gluing needs to preserve lengths of 



curves determined by two boundary value problems (see Definition 1.10). A technical point 



(which will be addressed in §4) is that one boundary component of one of the pieces (or 



more) in the decomposition will be singular, hence Theorem 0.4 cannot be applied directly. 

The paper is organized as follows. In §1 we introduce notation, recall a few useful facts 
concerning boundary value problems on graphs, and define a new natural metric induced 



by a solution of a boundary value problem (see Definition 1.9). In §2 we carry out analysis 



of level curves of the D-BVP solution; our study brings analysis and topology together in 
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order to provide a good notion of length for level curves as well as a topological description 
of them. In §3, we prove Theorem 0.4, and §4 is devoted to the proofs of Theorem |0. 1| and 



Theorem (12 Figures 2.10 - 2.12 , Figure 3.1 , and Figure 4.16 were generated by two lengthy 
programs written in Mathematica (version 7.0) by the author and will be available upon 
request. 



Remark 0.5. The assertions of Theorem 0.4 



may 



(in principle) be obtained by employing 
techniques introduced in the famous paper by Brooks, Smith, Stone and Tutte ([ID]), in which 
they study square tilings of rectangles. They define a correspondence between square tilings 
of rectangles and planar multigraphs endowed with two poles, a source, and a sink. They view 
the multigraph as a network of resistors in which current is flowing. In their correspondence, 
a vertex corresponds to a connected component of the union of the horizontal edges of the 
squares in the tiling; one edge appears between two such vertices for each square whose 
horizontal edges lie in the corresponding connected components. Their approach is based on 
Kirckhhoff 's circuit laws that are widely used in the field of electrical engineering. We found 



the sketch of the proof of Theorem 0.4 given in [TU] hard to follow. In fact, another proof of a 
slight generalization of this theorem was given by Benjamini and Schramm ([7], see also [5] for 
a related study) using techniques from probability and the dual graph of T. It is interesting 
to recall that it was Dehn [141 1903] who was the first to show a relation between square 
tiling and electrical networks. In an elegant paper, combining a mixture of geometry and 
probability, Kenyon ([23]) used the fact that a resistor network corresponds to a reversible 
Markov chain. He showed a correspondence between planar non-reversible Markov chains 
and trapezoid tilings. A completely different method, for the case of tiling a rectangle by 
squares, was given using extremal length arguments in [2E] by Schramm. One should also 
mention that Cannon, Floyd and Parry (see [12]), using extremal length arguments (similar 
to these in [26J), provide another proof for the existence of tiling by squares. Both [SB] 
and [12] are widely known as "a finite Riemann mapping theorem" and serve as the first 
step in Cannon's combinatorial Riemann mapping theorem ([11]). We include our proof of 



Theorem 0.4 , which is guided by similar principles to some of the ones mentioned above, 
yet significantly different in a few points, in order to make this paper self-contained. In 
addition, the important work of Bendito, Carmona and Encinas (see for example [4], [5], [6]) 
on boundary value problems on graphs allows us to use a unified framework to even more 
general problems. Their work is essential to our applications and we will use parts of it quite 
frequently in this paper as well as in its sequels ([19], [20]). 
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1. Preliminaries - boundary value problems on graphs 

We recall some known facts regarding harmonic functions and boundary value problems 
on networks. We use the notation of Section 2 in [3j. Let T = (V, F, c) be a finite network, 
that is a simple and finite connected graph with vertex set V and edge set E. Since in 
this paper Y = (V,E,c) is induced by (fl,dQ,T), we shall further assume that the graph 
is planar. Each edge (x,y) G E is assigned a conductance c(x,y) = c(y,x) > 0. Let 
V(V) denote the set of non-negative functions on V. If u G V(V), its support is given by 
S(u) = {x G V : u{x) 7^ 0}. Given F C V we denote by F c its complement in V . Set 
V(F) = {u G V{V) : S{u) C F}. The set dF = {(x,y) G E : x G F, y G F c } is called the 
edge boundary of F and the set <5F = {x G F c : (x, G E for some ?/ G F} is called the 
vertex boundary of F. Let F = F |J 5F and let F = {(x, y) G F : x G F}. Given F C V^, let 
r(F) = (F, F, c) be the network such that c is the restriction of c to F. We say that x ~ ?/ 
if (x,|/) G F. For x G F let A;(x) denote the degree of x (if x G 5(F) the neighbors of x are 
taken only from F). 

The following are discrete analogues of classical notions in continuous potential theory 

US]. 

Definition 1.1. ([4, Section 3]) Let u G V(F). Then for x G F, the function A«(x) = 
^ ri c(a;,j/) (u(x) — u{y)) is called the Laplacian of u at x, (if x G ^(F) the neighbors of x 
are taken only from F) and the number 

(1.2) E{u) = Y i &u{x)u{x) = c(x,y)(u(x)-u(y)) 2 , 

xeF (x,y)eE 

is called the Dirichlet energy of u. A function u G V(F) is called harmonic in F C V if 
An(x) = 0, for all x G F. 

For example when c(x, ?/) = 1, w is harmonic at a vertex x if and only if the value of u at 
x is the arithmetic average of the value of u on the neighbors of x. A fundamental property 
which we will often use in the maximum property, asserting that if u is harmonic on V C V, 
where V is a connected subset of vertices having a connected interior, then u attains its 
maximum and minimum on the boundary of V (see [27J Theorem 1.35]). 

For x G 5(F), let {yi,y2, ■ ■ ■ , y m } G F be its neighbors enumerated clockwise. The normal 
derivative (see [13]) of u at a point x G 5F with respect to a set F is 

(1.3) ^ (F)(X)= E c(x,y)Kx)- M (y)). 

The following proposition establishes a discrete version of the first classical Green identity. 
It plays a crucial role in the proof the main theorem in [12] and is essential in this paper. 

Theorem 1.4. ([3, Prop. 3.1]) (The first Green identity) Let F C V andu,v G V{F). 
Then we have that 

(1.5) ^ c(x,y)(u(x)-u(y))(v(x)-v(y)) = ^2Au(x)v(x)+ ^ g^( F )( x ) v ( x )- 

{x,y)£E x£F xeS(F) 



6 



SA'AR HERSONSKY 



Remark 1.6. In |3], a second Green identity is obtained. In this paper we will use only the 
one above. In [6] (see in particular Section 2 and Section 3), a systematic study of discrete 
calculus on n-dimensional (uniform) grids of Euclidean n-space is provided. Their definition 
of a tangent space may be adopted to our setting and does not require the notion of directed 
edges. 

Let T denote a fixed cellular decomposition of f2 U <9f2, and let F C be given. Let 
{c(x,y)}r xy \ e E be a fixed conductance function, and let T(F) be the associated network. 
We are interested in functions that solve a boundary value problem (D-BVP) on T(F). The 
following definition is based on [3], Section 3] and [5, Section 4]. 

Definition 1.7. Let k > be a constant. A solution of a Dirichlet boundary value problem 
defined on F(fi) is a function / £ "P(Q) such that / is harmonic in F, = and /|# 2 = 0, 
for some positive constant fc. 

Remark 1.8. The uniqueness and existence of a Dirichlet boundary value solution is guaran- 
teed by the foundational work in [3], Section 3] and [HI Section 4]. In fact, their work provides 
a detailed framework for a broader class of boundary value problems on finite networks. 

A metric on a finite network is a function p : V — » [0, oo). In particular, the length of 
a path is given by integrating p along the path (see [11] and [15] for a different definition). 
When p = 1, the familiar distance function on V x V is obtained by setting dist(A,B) = 
(J2xea 1) — 1 = A:, where a = (x, xi, . . . x^) is a path with the smallest possible number of 
vertices among all the paths connecting a vertex in A and a vertex in B. We now define a 
"metric" which will be used throughout this paper. 

Definition 1.9. Let F C V and let / e V(F). The flux-gradient metric is defined by 

(1.10) p{x) = ^(F)(x), if xG 5(F). 

This allows us to define a notion of length to any subset of the vertex boundary of F by 
declaring: 

(1.11) Length(5F) = | £ ?£(F)(x)|. 

In the applications of this paper, we will use the second part of the definition in order 
to define length of connected components of level curves of the D-BVP solution. In [THJ 
Definition 3.3], we defined a similar metric (/2-gradient metric) proving several length-energy 
inequalities. 
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2. TOPOLOGY AND GEOMETRY OF PIECEWISE LINEAR LEVEL CURVE 

Let Q be a polyhedral surface and consider a function / : — » R U {0} such that two 
adjacent vertices are given different values. Let v G G^°, and let wi, w^, ■ ■ ■ , be its k 
neighbors enumerated counterclockwise. Following (2U Section 3], consider the number of 
sign changes in the sequence {f(wi) - f(v), f(w 2 ) - f(v), . . . , f(w k ) - f(v), f(w{) - f(v)}, 
which we will denote by SgcJv). The index of v is defined as 



ffr) = 1 



(2.1) Ind 

Definition 2.2. A vertex whose index is different from zero will be called singular; otherwise 
the vertex is regular. A level set which contains at least one singular vertex will be called 
singular; otherwise the level set is regular. 

A connection between the combinatorics and the topology is provided by the following 
theorem, which may be considered as a discrete Hopf-Poincare Theorem. 

Theorem 2.3. ([2J Theorem 1], [2"l"l Theorem 2]) (An index formula) We have 
(2.4) J2 lnd f(v) = x(G)- 

Remark 2.5. Due to the topological invariance of x(G), note that once the equation above is 
proved for a triangulated polyhedron, it holds (keeping the same definitions for Sgc^(-) and 
Ind/(-)) for any cellular decomposition of x{G)- Also, while the theorem above is stated and 
proved for a closed polyhedral surface, it is easy to show that it holds in the case of a surface 
with boundary, where there are no singular vertices on the boundary (simply by doubling 
along the boundary). 

Henceforth, we will keep the notation of Section [0] and Section [TJ A key ingredient in our 
proofs of the theorems stated in the introduction is the ability to define a length for a level 
curve of g. The main difficulty in defining such a quantity (for level curves) is the fact that 
these are not piecewise linear curves of the initial cellular decomposition; hence we cannot 



directly define the weight p along them (see Equation (1.10)) 



Suppose that L is a fixed, simple, closed level curve and let 01,02 be the two distinct 
connected components of L in Q with L being the boundary of both (this follows by using 
the Jordan curve theorem; see for instance [HI Theorem III.5.G] for an interesting proof). 
We will call one of them, say Oi, an interior domain if all the vertices which belong to it 
have g- values that are smaller than the g- value of L. The other domain will be called the 
exterior domain. Note that, by the maximum principle, one of (9i,£?2 must have all of its 
vertices with g- values smaller than L. 

Let e G and x = e PI L. For x G" we have now created two new edges (x, v) and 
(u,x). We may assume that v G 0\ and u G Qi- We now define conductance constants 
c(v,x) and c(x, u) by 

(2.6) c(v, x) = < v ^)-9(u)) ^ = c(v,u)(9(u)-g(v)) ^ 

g(v)-g(x) g{u)-g{x) 
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We repeat the process above for any new vertex formed by the intersection of L with 
By adding all the new vertices and edges, as well as the piecewise arcs of L determined by 
the new vertices, we obtain two cellular decompositions, To 1 of 0\ and 7e> 2 of 2 . Also, two 
conductance functions are now defined on the one-skeleton of these cellular decompositions 



by modifying the conductance function for g according to Equation (2.6) (i.e. changes are 
occurring only on new edges). We will denote these by Co 1 and Co 2 respectively (on the arcs 
of L the conductance is identically zero). 

Definition 2.7. For i = 1,2, we let gi denote the solution of the D-BVP determined on the 
induced cellular decomposition on Oi defined by the following: 

(1) for any boundary component E of Oi and for every vertex x £ E, gi(x) = g(x), 

(2) the conductance function on Oi is Co,, and 

(3) gi is harmonic in O^ 



It follows from Equation (2.6), Definition and the existence and uniqueness theorems 



in [3] that for i = 1,2, gi exists and that gi = g\ot- We may now define a flux-gradient metric 
for L by using Equations (1.10) and (TTTTJ) . However, unlike the situation with the boundary 



components of fi, we have two possible choices, i.e. computing the normal derivative along 
L with respect to 0\ or with respect to 2 . Since in the applications we will cut Q along a 
particular L and wish to glue the resulting pieces together along L, these lengths computed 
with respect to the flux-gradient metric in each domain should be the same. The situation 
when L is not simple is more complicated and will be addressed after a detailed analysis 
of the topological structure of level curves has been carried out (see the discussion after 



Remark 2.19). This analysis is the main core of this section. The next lemma shows that 



all level curves are in fact closed. 

Lemma 2.8. A level curve for the function g is piecewise linear and closed, and each simple 
cycle of L contains at least one boundary component of dQ. 

Proof. The assertions of the lemma are certainly true for the components of dfl. Assume 
now that L is not one of the these level curves and furthermore that it is not closed. Let q 
be a boundary point of L. Since g is extended linearly on edges and in an affine fashion on 
triangles and quadrilaterals, q may (a priori) be an interior point of an edge of the cellular 
decomposition, a vertex, or in the interior of a cell. All of these cases are easily ruled out. 
It remains to prove that any such level curve contains at least one boundary component of 
dfl. If this is not the case, L (being a finite, closed, polygonal planar line) bounds a union 
of 2-cells in Q. This is a violation of the maximum principle for the harmonic function g\ 



constructed as in Definition 2.7| {0\ being one connected component of the union of the 



2-cells bounded by L). □ 

Remark 2.9. Since g is extended in an affine fashion along edges, it is clear that two disjoint 
level curves corresponding to the same g value are at combinatorial distance which is at least 
one. 

The structure of simple curves in the plane can be quite complicated and in fact is not 
completely understood. In our applications, we need only analyze the topological structure of 
closed curves arising as level curves of the affine extension (over triangles and quadrilaterals) 
of a harmonic function as defined above. We will henceforth work in the piecewise category. 
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Figure 2.12. The surface obtained 
by using g as a height function. 



Definition 2.13. A generalized piecewise bouquet of circles will denote a union of piecewise 
bouquets of circles where the intersection of any two bouquets is at most a vertex. Moreover, 
all such tangencies are exterior, i.e., no circle is contained in the interior of the bounded 
component of another. 

Proposition 2.14. Two simple cycles which correspond to the same g value either are 
disjoint or intersect at a single vertex. Furthermore, under these assumptions, a simple 
cycle is never contained in the annulus bounded by a different one (i.e. a "tangency" is 
always external). 

Proof. Let L\ and L 2 be two given simple cycles corresponding to the g- value I. Let A\ 
and A 2 be the two (piecewise) annuli which they bound respectively. Assume first that 
A\ PI A 2 = and that E = L\ fl L 2 is a piecewise arc. Let I and J be the endpoints of E. 



As mentioned before (see Lemma 2.8), both / and J are vertices. The link of the vertex / 
must contain a two cell (a triangle or a quadrilateral) such that an edge between two of its 
vertices crosses E. This implies that one of these vertices has g-value which is larger than 
I and the other has g-value which is smaller than I. This is absurd, since by the maximum 
principle all the vertices in A\ U A 2 which do not belong to g~ x {l) have g- values that are 
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strictly smaller than / (none of the vertices of these cells belong to g^ 1 ^)). It may also 
happen that L\ fl L 2 = {Pi, P 2 , • • • , Pc}, where k > 1 is an integer, and where the points Pi 
are oriented clockwise. Let t 2 be the arc connecting P\ to P2 which belongs to L 2 and which 
does not pass through any other of the Pi's. Let t\ be the arc in L\ which connects P\ to P2 
which does not pass through any of the Pj's. Then ti Ut 2 is a simple closed level curve which 



bounds an annulus (see Lemma 2.8). The link of the vertex Pi must contain a triangle or 
a quadrilateral with one of its edges crossing t\. As above, this leads to a violation of the 
maximum principle. 

Assume now that Aifl^ 7^ and (without loss of generality) that A\ULi contains an arc 
t of L 2 . Let its endpoints be P and Q (they both lie on L\). Then L\ U ((L2 \t) UP) consists 
of two simple cycles that satisfy all the conditions described in the case above; hence, this 
case may not occur. 

We now consider the case in which one of the annuli is contained in the other. Without 
loss of generality, assume that Ai C A 2 and that E = Li fl L 2 is a piecewise arc. Let / and 
J be the endpoints of E. The link of the vertex / must contain a two-cell, included in A 2 , 
where one of its edges connects a vertex in Ai to a vertex in A 2 \ (Ai U L{). This edge crosses 
L\] hence one of these vertices has g- value which is larger than I and the other has g- value 
which is smaller than /; this is absurd, since by the maximum principle all the vertices in A2 
which do not belong to g~ l {l) must have their g value strictly smaller than /. It may also 
be the case that L\ fl L 2 = {Pi, P2, • • • , P m } where m > 1 is an integer and the points p 
are oriented clockwise. Note that we allow all the Pj's to be a single vertex. Consideration 
of the link of any of the Pj's, and an argument similar to the one above, leads again to a 
violation of the maximum principle. 

Finally, assume that L\ flL 2 = 0. Observe that A 2 \A\ is an annulus; both of its boundary 



curves have g value I. By a construction analogous to the one described before Definition 2.7 
we obtain a new BVP problem defined on a cellular decomposition of (A2 \ A\) U L\ U L 2 . Its 
solution is harmonic on A2 \ A\ and has the same value on the boundary components; hence 
it is the constant function. The cellular decomposition of the annulus contains edges (part 
of edges) or vertices (and perhaps both) of the initial cellular decomposition. Moreover, the 
values of the solution and g coincide on (A2 \ A\) U L\ U L 2 . Both cases lead to the conclusion 
that connected pieces of edges of the initial cellular decomposition in the annulus have the 
same g-value, which is absurd. 

□ 

Theorem 2.15. Let L be a connected level curve for g. Then each connected component of 
L is a generalized piecewise bouquet of circles. 

Proof. First, we only use the fact that L is a closed polygonal line. Let S be the collection of 
all the self intersection points of L. By definition, edges of the polygonal (closed) curve L may 
cross only at points from S. If C is simple, the assertion follows immediately. Assume that 
L is not simple. It is well known that we may express L as a union of simple close polygonal 
curves. The intersection of any two cycles in this decomposition is a union of vertices and 



edges. Proposition 2.14 forces restrictions on any such decomposition. In particular, any 
two simple cycles in such a decomposition are either disjoint or intersect in a single vertex, 
none of which is contained in the annuli bounded by the other. The assertion of the theorem 
follows immediately. □ 
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Remark 2.19. Our interest here is in the existence of a decomposition of the level curve into 
simple pieces. One may also consider the question, or rather a series of questions, regarding 
a decomposition of the non-simple polygon determined by a general closed polygonal curve. 
For many deep questions and an excellent survey on this interesting subject (some of the 
questions involved are known to be NP), see for instance [22] . 

The topological structure of a general connected component of a level curve L (Theo- 
rem 2.15) means that the complement of L in Q U dQ is composed of a bounded domain 
(which is a disjoint union of annuli) in which all vertices have ^-values smaller than the g- 
value of L, and a second domain whose boundary consists of L and possibly other curves in 
dQ), where L is in the boundary of both. Following a construction similar to the one preced- 



ing Definition |2.7[ we may now define two modified conductance functions and two solutions 
of D-BVP problems. By abuse of notation, we call the first domain 0\, the second O2, 
and we let C^, Co 2 , gi, #2 be the analogous quantities. The corresponding two new cellular 
decompositions and their union will henceforth be called the induced cellular decomposition. 

Theorem 2.20. Let L be a level curve for g. Then the following equality holds 

(2.21) l 0l (L) = l 02 (L), 

where the left-hand side denotes the length of L measured with respect to the flux-gradient 
metric induced by g\, and the right-hand side denotes the length of L measured with respect 
to the flux-gradient metric induced by g 2 . 
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Proof. First, let us assume that L is simple and corresponds to the g- value h. Since L is at 
combinatorial distance at least one from any other level curves of g of the same value, the 
1-combinatorial neighborhood of L in 2 is comprised of vertices that are of ^-values greater 
than h. 

Let x G fl L. Let yi, y 2 be its neighbors in L, {x\, . . . , x m } in 0\ and {zi, . . . , in 
(!?2- Since g is harmonic at x we have 

(2.22) = Ag(x) = £ c(x, y) (</(a;) - </(j/)) . 
Hence, (since flr(ar) = ^(^/i) = ^(2/2)) 

m k 

(2.23) = c ( x > x i) - g{xi)) + ^2 c(x, zj) (g(x) - g(z j )) . 

i=i 3=1 

Let x G 7c?! U 7o 2 be a new vertex in L (that is x ^ T {0) ) with u G Ci and u G C 2 - We 
have 



(2.24) 

and 

(2.25) 



dg 



dn 



dg 2 
dn 



-{O x ){x) = c(v, x)(g(x) - g(v)) = c(u, v)(g(u) - g(v)), 



(0 2 )(x) = c(u,x)(g(x) -g{u)) = c(u,v)(g(v) - g(u)). 



By summing the three equations above over all vertices in L, the assertion follows. 

By combining the arguments above with the topological structure of a general level curve 
provided by Theorem 2.15 and a simple modification at singular vertices along L, we will 
now treat the second case. Let us assume that L is not simple. Recall that a level curve 
may intersect itself only at a vertex. Let S = {fi,. . . ,v m } be all such intersection points 
and {Ind/(ui), . . . , Ind/(t> m )} be their indices respectively. For a vertex x G L (of the initial 
triangulation or the induced one) which is not in S (note that Ind/(x) = 0), the arguments 
above leading to Equation 2.22 go through precisely in the same manner and yield (with 
Ox, O2 modified as explained immediately after Remark 2.19) the same conclusion. 

Let v G S. Then v is in the intersection of C±, . . . ,Ci piecewise simple circles where, / = 
Sgc g (v )/2 circles. Let {x\, . . . , x\^} be the neighbors of v in 0\ fl interior(Cj), . . . , y m } 
the neighbors of v in L and {z\, . . . , z p } in 2 - Since g is harmonic at v , we may now modify 
slightly Equation (2.23). Since g(v) = g(yi) = . . . = g(y m ), we obtain 



(2.26) 



1 

= J2J2 C ^ X ^ fo( 

3=1 i=i 



v 



g[x; 



')) + '52 c ( v > z j)(9{v)-g{zj)). 



The assertion of the theorem now follows by summing over all the vertices in L. Finally, if L 
has several connected components, the assertion of the theorem follows by addition over each 
component (which must be at combinatorial distance greater than one from each other). 

□ 



Remark 2.27. We now wish to start using Theorem 2^ and Remark 2J5 One should note 
that since boundary components are at the same g-level, a modification is needed. This is 
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done in the following way. Along each boundary component, add a piecewise linear curve 
homotopic to it, and at combinatorial distance which is smaller than one. This can be done 
in such a way that there is exactly one sign change at each new vertex (not taking into 
account vertices on the boundary). Let Q G Q be the region bounded by these curves. 
Doubling along the boundary of Cl gives a surface without boundary with all of its singular 
vertices identical to those of g, for which the above mentioned theorem may be applied. 

We end this section with a useful property of level curves. This property is essential in 



the proofs of Theorem 0.1 and Theorem 0.2 In the course of these proofs, we will need to 
know that there is a singular level curve which encloses all of the interior components of 
dfl. This level curve (which is necessarily singular) is the one along which we will cut the 
surface. In the inductive process, we keep cutting along such level curves over domains of 
fewer boundary components until the remaining pieces are annuli. 

Proposition 2.28. With notation as above, there exists a singular level curve which contains 
in the interior of the domain it bounds, all of the inner boundary components of dVt. Such 
a level curve is unique. 

Proof. We argue by contradiction. Throughout the proof, we use the topological structure of 



level curves provided by Theorem |2.15[ Suppose that the first assertion of this proposition 
is false. Then every singular level curve omits at least one inner boundary component. Let 
L be such a curve. Suppose that L omits k > 1 boundary components and contains I > 2 
boundary components in its interior, where k + l = m—l. Let {v i, v 2 , ■ ■ ■ ,v n } be the singular 
vertices that belong to L; let {wi, W2, ■ ■ ■ ,w m } be the singular vertices in the interior of L. 
It is easy to see that (for example by induction) 

(2.29) Yl lnd g( x ) = 1 - l 

xe{uoi,...,v n } 



By Theorem 2.3 we must have additional singular vertices {qi, . . . ,q s } such that 

(2.30) lnd M = (!-*)-!■ 

ye{qi,...,q s } 

Let L\ be a singular level curve that has q\ as a singular vertex. Suppose that L\ does not 
enclose L in its interior. Then 

(2.31) Ind 9 ( gi ) < (1-k), 

with equality if and only if L\ encloses all k boundary components. Therefore, we must have 
s = 2, and a singular level curve that passes through q 2 will necessarily enclose all interior 
components of dQ. Hence, we arrive at a contradiction. 

Assume now that L\ encloses L and t > 1 out of the k boundary components. Let 
{qi, q 2 , . . . , q s } be the singular vertices on the various cycles on L\ which do not belong to 
the cycle containing L or its interior. It is easy to see that 

(2.32) Ind s (z) = (l-*). 

z&{qi,q2,---,q s } 



By Theorem 2.3 we must have additional singular vertices {pi, . . . ,p n } such that 
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(2.33) Yl Ind ff (*) = -l-(Jfe-t). 

te{pi,...,q„} 

We now repeat the process above with Li replacing L, l + t replacing I and k — t replacing k. 
After finitely many times (at most k — 1), we must end up with a singular level curve with 
a single singular vertex of index —2 that encloses all inner boundary components of dQ. 

□ 



We are now concerned with geometric properties of the set of lengths of level curves. This 
is significant for the applications. For example, even in the case of an annulus (which is the 



subject of Theorem 0.4), it is not clear that level curves in the same homotopy class have the 
same length. Since we wish to map a given annulus to a straight cylinder, we must verify 
this property and its suitable generalizations. 

Theorem 2.34. There exists a finite set of non-negative numbers K, such that the following 
holds: 

(1) /C contains ; and k is its maximal element, 

(2) /C is monotone decreasing, 

(3) any level curve corresponding to a g-value which does not belong to K is simple, and 

(4) any component of a level curve, which corresponds to a g-value which is strictly 
between any two values k n > k m in IC, such that k n — k m is minimal, and is contained 
in a unique simple cycle C n determined by k n , has its l\-length equal to that of C n . 

(5) Moreover, the length of C n is equal to the length of the component of g~ 1 (k m ) which 
it encloses. 

Proof. We let K, denote the set of critical values of g union k and 0; this gives assertion 
(1). Since x(^ U 9Q) is negative and the index of each singular vertex is less than or equal 



to —1, Theorem 2.3 implies that the number of singular vertices, and hence the number of 



corresponding critical g-values, is finite. Assertion (2) follows by ordering. Assertion (3) is 



the content of Lemma 2.8 We now turn to proving (4). Let L be a connected component 



of a level curve as described in (4). By the structure theory provided in Theorem 2.15| and 
the maximum principle, L is contained in a unique simple cycle, part of the bouquet which 
composes g~ 1 (k m ). Let C n denote this cycle. We first assume that the combinatorial distance 
between L and C n is greater than or equal to 2. We apply the construction that precedes 



Definition 2.7 to the level curves C n and L and the annulus Al,c„ enclosed by them. By 
abuse of notation, we will keep denoting the modified set of vertices by . 

Let F = Fc n ,L denote the subset of fl Al,cv Since the combinatorial distance between 
C n and L is greater than or equal to 2, it follows that 



(2.35) F ^ and that 5{F) = {v G T (0) : v G C n U L}. 

Recall that g\, the modified solution of the D-BVP, is defined on F U 5(F) by requiring 
that gi is harmonic (with respect to the modified conductance function) in F, gi\i = g\h 
and gi\c n — g\c n - By the uniqueness of the D-BVP solution, it is clear that g\ = g\p. Let 
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F° denote the set of vertices in F which do not belong to 5(F). We apply Proposition 1.4 
(the first Green Identity) to F with u = g\ and v = 1. Equation (1.5) then gives 



(2.36) 

Hence 
(2.37) 

which implies that 
(2.38) 







xe<5(F) 



dn 



E?>°>M+Et^°)< 



a:6i 



It follows that, with respect to the flux-gradient metric (see Equation (1.11 )), L and C n have 



the same length. We now deal with the case that the combinatorial distance between C n and 
L is less than 2. While we wish to use again the first Green Identity, care is needed since F 
is empty in this case. We add a vertex on each edge in the modified cellular decomposition 
that has a vertex v on L and a vertex u on C n . The value of g\(t) is defined by the value of 
gi on this edge and we also let 



(2.39) d(v,t) 



c(v,t)(g 1 (v) -gi(u)) 

gi(v) -gx(t) 



and c(u, t) 



c(u,t)(g 1 (u) - gi(v)) 
gi(u)-gx(t) 



where c(v,t) and c(u,t) are defined as in equation (2.6) with t replacing x. We keep all 
conductance functions on other edges unchanged. 

By applying again the existence and uniqueness theorems in [3], we obtain a solution g 2 
of a new D-BVP defined on Al,c„ by requiring that g 2 \i = 9i\l, 92\c n = 9i\c n an d that g 2 
is harmonic in Al,c h - ^ follows that g 2 is g% on all vertices in and is modified so as to 
have the values of g\ on all vertices defined above. It is easy to verify (using Equation 2.6 
and Equation 2.39) that 



(2.40) 



()92 (A L , Cn )(v) = i£(A L , Cn )(v) and d -B(A L>Cn )(u) = ^(A L , Cn )(u). 



on on on on 

Assertion (5) is proved by following the same techniques that were used in proving assertion 



(4) and Theorem 2.15 



□ 

Remark 2.41. Observe that once assertion (5) is proved, in order to prove (4), one needs only 
add at most one vertex on each edge whose vertices lie on C n and g~ 1 (k m ) respectively. Also, 
the same techniques employed in the proof of assertion (4) allow one to get equality between 
the sum of the Zi-lengths of all the components for a non-singular level curve enclosed in a 
simple cycle and the Zi-length of this cycle. Finally, the case g~ x (0) was not stated separately, 
yet it similarly gives that the sum of the /i-lengths of the inner boundary components of 
dQ equals the /i-length of the outer boundary. We will revisit this last point in the proof of 
Theorem 10.11 
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3. THE CASE OF AN ANNULUS 



In this section we prove Theorem |0.4| The proof consists of two parts. First, we will show 
that there is a well-defined mapping from into a set of (Euclidean) rectangles embedded 
in the cylinder S^. The crux of this part is the fact that level curves for g have the same 



induced length (measured with the l\ metric, see Remark 2.41 ) and a simple application of 



the maximum principle. Second, we will show that the collection of these rectangles forms 
a tiling of with no gaps and no overlaps. This is a consequence of the first part and an 
energy-area computation. The latter follows by our construction, the dimensions of S_a, and 



the first Green identity (see Theorem 1.4). 




Figure 3.1. Level curves 
in an annulus 

Given a straight Euclidean cylinder of height h, we will endow it with coordinates arising 
naturally from S 1 x [0,h]. The boundary component corresponding to h will be called the 
top and the other one will be called the bottom. Before providing the proof, we need a 
definition which will simplify keeping track of the mapping /. 

Definition 3.2. A marker on a straight Euclidean cylinder is a vertical closed interval which 
is the isometric image of 9 x [a, b], for some 9 G [0, 2tt) and [a, b] C [0, h] with a < b. The 
marker's uppermost end-point corresponds to {9,b) and its lowest end-point to (9, a). 



Proof of Theorem 0.4, Let be a straight Euclidean cylinder with height H = k and 
circumference 



C 



dn 



Let £ = {Li, . . . , Lk} be the level sets for g corresponding to the vertices in arranged 
in a descending g-values order. We place a vertex at each intersection of an edge with an 
Li, i = 1, . . . , k,, and if necessary more vertices on edges so that any two successive level 
curves in £ are at combinatorial distance (at least) two. We call the first group of added 
vertices type I and the second type II (recall that conductance along edges are changed as 



well, according to the discussion preceding Definition 2.7). In particular, the assertions of 



Theorem 2.34 and Remark 2.41 hold. Since x(*A) = and the index of a singular vertex is 
always negative, the index formula (Theorem 2.3[ ) prevents the existence of singular vertices 
in A. Therefore, /C = {0, k}, and furthermore the length of any g- level curve is equal to C. 
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Let the vertices {x±, . . . ,x p } in E\ = L\ be ordered counterclockwise (and on any level 
curve), starting with x\. Let {yi,y2, ■ ■ ■ ,Vt} be its type I neighbors in the induced cellular 
decomposition oriented clockwise (which will always be the ordering for neighbors). We 
identify x\ with x k in the coordinates above. We associate markers {m Xuyi , . . . ,m Xuyt } 
with X\ in the following way. The length of marker m xit y s , for s — 1, . . . , t, is equal to (the 
constant) g(xi) — g(y s ) and its uppermost end-point is positioned at the top of at 

(3.3) Xl + T.TM^m)(9{^) - g{y k )) x 

measured counterclockwise. For each edge e UjV = [u,v] with g(u) > g(v), let Q u>v be a 
Euclidean rectangle with height equals to g(u)—g(v) and width equals to c(u, v)(g(u) —g(v)). 
We will identify a planar straight Euclidean rectangle and its image, under an isometry, in 
a straight Euclidean cylinder. For s — 1, . . . ,t, we position Q xi ,y a in Sa m sucn a wa Y that 
its leftmost edge coincides with m xuy3 . By construction and the position of the markers, 



Qxi,y s l~l Qxi,y s+ i m xi,y s+ i- 

Assume that we have placed markers and rectangles associated to all the vertices up to x k 
where k < p; let Z\ be the leftmost neighbor of x^+i and let Q Xk ,v be the rightmost rectangle 
associated with x k - We now position the marker m Xk Z1 , associated with Xk+\ and z\ so 
that it is lined with the rightmost edge of Q Xk ,v and his upper end-point is at the top of S^. 
We continue placing markers and rectangles corresponding to the rest of the neighbors of 
Xk+i- We terminate these steps when k = p. Note that the top of Sa is completely covered 
by the top of the rectangles constructed above where intersections between any of these (top 
edges) is either a vertex or empty. 

For all 1 < n < k, assume that all the markers corresponding to vertices in L n _i and their 
associated rectangles have been placed as above in such a way that the following condition, 
which we call consistent, holds. For [w,v] G with g(w) > g(v) and s e [w,v] a vertex of 
type I, the uppermost end-point of the marker m StV coincides with the lowest end-point of 
the marker m WiS ; moreover, the two rectangles Q WyS , Q s , v tile Q VjW . Informally, this condition 
allows us to "continuously extend" rectangles associated with edges that cross level curves 
along these, and therefore will show that that edges in are mapped (perhaps in several 
steps) onto a unique rectangle. 

We will now show how to place the markers and rectangles corresponding to the vertices 
of the level set L n in a consistent way. The uppermost end-point of each marker associated 
with a vertex v G L n in this level set, and any of its neighbors in L n+X is placed in S^ 
at height g{v). Observe that v is a vertex in some [qi,v], where belongs to L n . Choose 
among all such edges the rightmost (viewed from L n ). Let [qo,v] be this edge and let m q0tV 
be its marker. Place the marker of v which corresponds to an edge [v, w] with w being the 
leftmost vertex among the neighbors of v in L n+ ±, so that its uppermost end-point coincides 
with the lowest end-point of m q0tV . Let s be a vertex of type I on L n . By definition, sq 
is connected to a unique vertex vq G L n+ i and to a unique vertex wq G L n -\. Let s q be 
the first vertex to the left of s of type I in L n . By definition, s q is connected to a unique 
vertex w p in L n _i and to a unique vertex vi G L n+1 . Let {si, . . . , s 9 _i} be the vertices on L n 



(3.4) 
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between so and s q and let {wi, . . . , w p -\} be the vertices on L„_i between u> and u> p . Let 

Ql = Q^O, 

SQ,s q ,w p be the piecewise linear rectangle enclosed by [iwoj ^o] U [wpj ^g] U ^n—x U -^n? 
and which contains {so, . . . , s 9 }, and let Q 2 = Q So ,v ,s q ,vi be the piecewise linear rectangle 
enclosed by [sq, t>o] U [s 9 , uj] U L n U L n , and which contains {so, . . . , s q }. 

In order to prove that the consistent condition holds for all markers and rectangles created 
in this step, it suffices to prove it at s q , assuming (with out loss of generality) that the 
first marker we placed in a consistent way is m SO)Vo . By construction (see in particular 



Equation (3.3)) we need to prove that 
(3.5) 



9-1 



dg 

+ ^(22) ri ght(so) 



where the subscript "right" indicates that neighbors in the expressions are taken from Qi 
(first line) or Q 2 (second line) only. It is easy to check that since the modification of g is 
harmonic at each Si,i = 0, . . . q, and since so and s q are type I vertices, Equation (3.5) holds. 

To conclude the construction, continue as above, exhausting all all vertices in L n . By 
the maximum principle, our construction, and the fact that all level curves have their flux- 
gradient length equal to C, it is clear that the union of the rectangles is contained in S^. 


















4 



Figure 3.6. Several markers and rectangles on the Euclidean cylinder after 
the completion of the construction. 



We now prove that the collection of rectangles constructed above tiles S4 leaving no gaps. 
Without loss of generality, suppose that the collection of rectangles does not cover a strip 
of the form [9 1 ,9 2 ] x [h ,hi] in SU, where < 9 1 < 6 2 < 2n and < h < hi < k. By 
harmonicity, there exists at least one path whose vertices belong to such that the values 
of g along this path are strictly decreasing from k to 0. 

In particular, the value h\ is attained on some edge or vertex of this path. By construction, 
a gap in a (7-level curve (i.e. an arc of level curve which is not covered by the top edges of 
rectangles) never occurs when h\ is the (/-value associated to a vertex in the modified cellular 
decomposition. Hence, we may assume that h\ is attained in the interior of an edge. Let 
Lf n be the corresponding level curve. Recall that L/ ll is simple and closed (as are all other 



level curves of g\ see Lemma 2.8). We now follow the construction preceding Definition 2.7 
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and let {ui,Uz, . . . ,u q } be all the new vertices on L^, that is we place a vertex at each 
intersection of an edge in with L^. As mentioned in the first paragraph of the proof, 
the flux-gradient length of is equal to C. Moreover, this length is equal to 

(3-7) E§^)(^)> 

?.=i 

where 0\ is the interior of the annulus enclosed by and E% (see the discussion proceeding 



Definition 2.7). In particular, we may now place markers and rectangles associated to the 
collection of vertices {u±, u 2 , ■ ■ ■ , u g } so that L hl is covered by the top edges of these rectan- 
gles. Recall that up to this moment, markers and rectangles were associated only to vertices 
of type I in addition to those in T^°\ It is very easy to check that this can be done in a 



consistent way (using Equation (2.6)). Since g is extended affinely over edges, every value 
between ho and h\ is attained by g. Repeating this argument shows that all level curves are 
covered by rectangles. Hence the collection of rectangles leaves no gaps in S^. 

Using an area argument, we now finish the proof by showing that there is no overlap 
between any two of the rectangles. Let TZ be the union of all the rectangles. By definition, 

(3.8) Area(ft) = c(x, y)(g(x) - g(y))(g(x) - g(y)). 



Note that the sum appearing in the right-hand side of Equation (3.8) is computed over 
the induced cellular decomposition. A simple computation (using again Equation (2.6)) 
and the fact that the construction is consistent shows that this sum is equal to the one 
taken over [x,y] E Hence, the right-hand side of this equation is the energy E(g) 



see Definition 1.1[ ). Therefore, by the first Green identity, applied with u = v = g (see 



Theorem 1.4), and the dimensions of S^, we have 
(3.9) E(g) = Area(ft) = Area(SU). 

Hence, since the rectangles do not overlap, they must tile S4. It is also evident that the 
mapping / constructed above is energy preserving. 

□ 

4. THE CASE OF AN m,- CONNECTED BOUNDED PLANAR REGION, m > 2 



There is a technical issue which we need to address in the proofs of Theorem |0.2| and The 
orem 



0.1 Since Proposition 1.4 is applied all along, we must make sure for example, that the 
combinatorial distance between any two distinct level curves is at least two. Therefore, one 
adds as many vertices of type I and type II as needed, and modifiy conductances accordingly 



(as done at the beginning of the proof of Theorem 0.4) until the assertions of Theorem 2.34 
and Remark 12.411 hold. 



4.1. Pair of pants. In this subsection we provide a proof of Theorem |0.1| One interesting 
ingredient of the proof is the construction of a good mapping (in the sense of Theorem 0.4) 
from a planar annulus with one singular boundary component into a Euclidean cylinder 
having one boundary component, a singular curve and the other one which is simple. The 
construction of the mapping / is then achieved by gluing two Euclidean cylinders to the one 
above along the singular component. 
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Proof of Theorem 10. 11 

We extend the solution g to an affine map g : \V U dP\ — > R + U and use g as a height 
function for \V U dP\. This gives a two-dimensional polyhedron, denoted by P g , which is 
homotopically equivalent to V. We have 

(4-1) x(P g ) = X(V) = -1. 

Double P g along its boundary to obtain a genus 2 closed surface which will be denoted by 
DP. The Index Theorem (see Theorem |2.3[ ) asserts that 



(4.2) 



E 

veDv 



lnd g {v) 



-2. 



It follows from Equation (2.1) and a simple Euler characteristic computation that we must 
have a single vertex u G P g whose index equals —1, which means that Sgc(w) = 4. Let L(u) 
be the figure eight curve corresponding to the value g(u). Then L(u) = Lj U L n , where Lj 



and Lji are two simple cycles intersecting at u (using here the assertion of Theorem 2.15). 

Moreover, since E\ and E\ are the only inner components of <9f2, by applying Lemma |2!8| 
we may assume that E\ is contained in Li and E\ is contained in Ljj. By cutting along L(u), 
we can decompose Q U dQ into three components. The first, denoted by Q\, has boundary 
components E\ and L(u). The second, denoted by f^, has boundary components E\ and L/, 
and the third, denoted by Q|, has boundary components E\ and Ljj. Note that the interiors 
of these components are all homeomorphic to S 1 x (0, 1), and moreover that Ql, and Vl\ are 
each homeomorphic to an annulus, however — — 1- Therefore, Qi may be viewed as 

an annulus with one singular boundary component, or equivalently as a really short pair of 
pants, i.e. one with no cuffs. We now define three D-BVPs on the above domains following 
the discussion preceding Definition 2.7 Note that an essential property of L(u), which is 



given by Theorem 2.15 allows us to define a D-BVP on fix- 



Let gi be the solution on f2i, g\ on and g\ on Vt\. By Remark 2.41, we have that the 



lengths of Lj and E\, with respect to the flux-gradient metric induced by g\, are both equal 
to 

(4.3) ^ d9 



Similarly, we have for that the lengths, with respect to the flux-gradient metric induced 
by g\ of Lji and £f , are equal to 

(4.4) <^ d9 



Finally, by applying Theorem 2.34 with C n = E\ and L = L(u), we have that the lengths 
of Ei and L(u), with respect to the l\ gradient metric induced by g±, are equal to 

(4.5) 



E %m 



We now turn to the construction of (S-p, f). First, apply Theorem 0.4 to ^2 and to Q 



The outputs are two straight Euclidean cylinders S n i and S n 2, two mappings f n i : fi^ — > <S\ 
and /fj2 : Vt\ — > S2 (where the maps, their domains and the cylinders are described in detail 
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E 1 E 2 

2 2 

Figure 4.6. The construction of the Euclidean pair of pants with one singular point. 



in Theorem 0.4). Second, we wish to apply Theorem 0.4 to However, one boundary 



component of it (L(u)) is singular at u. Note that the construction described in the proof of 



Theorem 0.4 will go through for any annulus contained in Qi with one boundary component 



being E\ and the other being L g (L(u))+e (with e > arbitrarily small). Thus, we obtain a 
sequence of Euclidean cylinders that have the same circumference (the flux-gradient metric 
length of Ei) and their heights converge to k — g(L(u)). 

As e — > 0, the level curves corresponding to g(L{u)) + e converge geometrically to L{u). 
Hence, the limiting cylinder Sq is a Euclidean cylinder with one singular boundary compo- 
nent (bottom). We call one boundary component of this cylinder top and the other bottom. 
Concretely, Sq is obtained by taking a Euclidean cylinder of height k — g(L(u)) and circum- 
ference which is equal to the flux-gradient metric length of E\, picking two points on one 
boundary component of it at a distance which equals the flux-gradient metric length of Li, 
and identifying them. We will abuse notation and denote this point by u. Since identification 
occurs only at the boundary, the tiling persists in the limit. We denote by /q 1 the mapping 
given by Theorem 0.4[ modified on the bottom as above. Let 



(4.7) 



/ = U / n i U /, 



where the domain of / is U U ^2 an d its image is obtained by gluing isometrically 
Si and S 2 to the bottom in such a way that the intersection consists of only one point (see 
Figure 4.6), denoted by f(u). By Remark 2.41, we have that the length (measured with 



the Zi-induced metric) of the bottom of Sq is equal to the sum of the lengths of the tops 



of the two cylinders Si and S2 (one can also obtain this by applying Proposition 1.4 and 
Equations (4.3 )-( 4~5| ) directly). Furthermore, we require the gluing to be consistent; i.e the 
gluing described above respects rectangles corresponding to edges which cross L(u) (see the 



Proof of Theorem 0.4). The fact that this can be done may be justified by basically the same 
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argument appearing before and after Equation (3.5). Note that one consequence of this is 



that f(u) is uniquely determined (the modified g being harmonic being the key issue) 
It is easy to check that the only cone angle is at this point and is equal to An. 
mappings / n i 



fn 2 2 



0.4 



The 

the mapping 



are energy preserving in the sense explained in Theorem 
/fij is also energy preserving (since the identification is done on bottom only) . Since the 
gluing is by isometries, and the cylinders involved meet only at one vertex, the mapping / 
is energy preserving. 

□ 



4.2. The case of an m-connected domain, m > 3. 



of Theorem 0.2 The proof is a modification of Theorem 0.1 



In this subsection we provide a proof 
with some bookkeeping and 



successive changes of the initial D-BVP. 

Proof of Theorem 0.4. Let L\ be the unique level curve whose existence is guaranteed 
by Proposition 2.28 Let Q\ be the domain bounded by E\ and L\. Observe that f2i is 
homeomorphic to S 1 x (0, 1). Suppose that L\ = Zi U 1% . . . U Z&, where each Z, is a simple 
cycle (for % = 1, . . . ,k), where any t wo of these cycles are either disjoint or intersect at a 
singular vertex in (see Theorem 2.15). Let U = {ui, . . . ,u m } be the set of all singular 
vertices along L\ with {Ind s (wi), . . . , Ind 9 (w m )} being their indices, respectively. 

The interior of each Zi, which will be denoted by Z°, is a gi-connected domain with qi < 



ni 



1, for 



1, . . . , k, (unless k = 1). We now modify the initial D-BVP on Q\, U and Z° 
as described in the discussion preceding Definition |2.7[ and we obtain k + 1 new harmonic 
functions: go which is defined on E\, Qi and L\, and gx,...,g^ defined on li, . . . , Z& and their 
interiors, respectively. 

Recall that modifications are made by adding vertices to along the Zj's and by defining 
conductance along new edges: those which cross the Zj's (See Equation (2.6)). By part (5) 
of Theorem 2.34 and Remark 2.41 we have that the length of Li, with respect to the flux- 
gradient metric induced by go, is equal to the length of E\ with respect to the flux- gradient 
metric induced by go. Hence both are equal to the length of E\ with respect to the initial g 
metric. We record this as 

(4.8) lengthy (E0 = lengthy [L x ). 

By using in the interior of each Zj, we now choose a singular curve enclosing all of its 
5i — 1 boundary curves (see Proposition 2.28). The result is a set of singular level curves 
W = {Wi, . . . , W n } and the set of singular vertices V = {v i, . . . , v p } they contain with their 
associated indices X = {Ind g (vi), . . . , Ind 9 (w p )}. Each simple cycle of a singular level curve 
Wj contains Sj < qi boundary curves unless W° is an annulus. As above, we have only added 
vertices on the Zj's, assigning each one of the vertices on a specific Zj, the g(Zj)-value. The 
conductance constants are changed only along new edges, those which cross Zj, according to 



Equation (2.6). Hence, the set W coincides with the set of singular level curves of g minus 



Li (similar statements hold for V and X). 

We repeat this process (at most) finitely many times, modifying (if needed) the cellular 
decomposition and defining conductance constants according to Equation ( |2.6 ), until the 
interior of each simple cycle of each singular level curve is an annulus. At each step, we 
obtain new harmonic functions defined on domains with fewer boundary components. An 



equation analogous to Equation (4.8) holds for each simple cycle which is a component of 
a singular level curve and the nearest singular level curve it contains. That is, its length 
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measured by the flux-gradient metric induced by the harmonic functions defined on its 
interior is equal to the length of the singular level curve, both measured by the harmonic 



map defined in the interior of the simple cycle (see part (5) of Theorem 2.34). In turn, they 
are equal to the length of the simple cycle measured by the harmonic function defined on its 
exterior (see Theorem 2.20). 

It is evident from the proof of Proposition 2.28 the disjointness of level curves which 
correspond to different g-values, and the maximum principle that there is a well-defined 
hierarchy on the set of singular level curves. Each component of a singular level curve, say 
g l {k m ), k m e /C, is contained in the interior of a domain whose one component is a unique 
simple cycle. This cycle belongs to a unique singular curve, say g _1 (fc n ), k n e /C, k n > k m (or 
Ei). Simple cycles corresponding to level curves are either contained in a domain bounded 
by a simple cycle which is part of a singular level curve, or are parts of singular level curves 
as is detailed in Theorem 12.151 

We now turn to the construction of Sq and the mapping /. We start with a straight 
Euclidean cylinder, Ce 1 ,l 1 i of height k — g(Li) and circumference which is equal to 



(4.9) 



dn 



Since L\ is a generalized bouquet of circles, we can select a finite number of points in the 
bottom and identify subsets of them in such a way that the quotient is topologically and 
metrically isomorphic to L\. 




Figure 4.10. A quotient of a boundary component of a Euclidean cylinder 
yielding a generalized bouquet of circles 



As in the proof of Theorem 0.1 we wish to apply Theorem |0. 4 to Vt\^dVt\ {dVt\ = ExULi). 
Since L\ is singular, we may not apply it directly. Instead, we follow the argument carried out 
in the proof of Theorem 0.1 All level curves (which are simple) corresponding to g{L\) + e, 
with e > being small, will have the same length as L\ (see part (5) of Theorem 2.34); hence, 
as e — > 0, the sequence of Euclidean cylinders (guaranteed by Theorem 0.4 applied to the 



sequence of annuli bounded by these level curves and E{) with their tops being fixed, heights 
equal to k — (g{Li) + e), and circumferences equal to YlxeE ^{ x )i converges to Ce 1 ,l 1 with 
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the top persisting in the limit and the bottom replaced by its quotient as explained above. 
We denote the mapping given by Theorem 0.4 and modified as above on the bottom by fn 1 . 



For each U, i = 1, . . . k, there are two cases to consider. First, suppose that 1° does not 
contain any singular level curve. Without loss of generality, let the unique boundary curve it 



contains be E\. By assertion (4) of Theorem 2.34, we conclude that the lengths of E\ and I 



are equal (measured with the ^-induced flux-gradient metric). Thus, Ai = 1° U U U E\ is an 



annulus with two (non-singular) boundary components. We may therefore apply Theorem 0.4 
and obtain a mapping fi '■ Ai —> C lijE % (= S^J. We now attach the top of the resulting 
Euclidean cylinder Ci uE i to U by an isometry which is consistent. We conclude that the 
gluing may be taken as such, by first arguing that the length of U measured with respect 
to the flux-gradient metric induced by g , is the same as its length with respect to the 
length induced by the flux-gradient metric of <ft. This is justified directly by employing 



Theorem 2.20 or by using a series of equalities similar to those in Equations (4.3)-(4.5) and 
the reasoning leading to them. The fact that the gluing is consistent follows from the same 
arguments we have used before. 

Second, assume that 1° contains a singular curve, say Wi G W. Recall that Wi is the 
unique singular curve in 1° which encloses all of the boundary curves in 1°. With U replacing 
Ei, Wi replacing L\, and gi replacing g , we repeat the construction, yielding a Euclidean 
cylinder with its bottom being a singular curve as described above (in the case of Li), and a 
mapping which we will denote by fi it w v The cylinder Ci it Wi * s attached to Ce 1 ,l± by gluing 
its top to the simple cycle corresponding to U in the quotient of the bottom of Ce u Li, that 
is, the singular component of Ce 1 ,l 1 - 

The above analysis is now carried out, at most finitely many times, for each li, i = 1, . . . , k, 
until we are left with Euclidean cylinders having both of their boundary components non- 
singular. In particular, the bottom of each cylinder is the image (under the map given by 



Theorem 0.4) of E l 2 and its length is equal to 

dg_ 
dn 



(4.11) - f^O*)' for * = V" >™- L 



By construction, the images of E\ and £|, i — 1, • • • , m, under the maps whose construction 
is described above, are the only boundary components of Sq. The map / is the obvious 
union of the collection of mappings constructed above and is analogous to the one appearing 



in Equation (4.7). It is also evident that / is energy preserving in the sense described before. 

We now compute the cone angle, 4>{y), at a singular vertex v with Ind g (t>) = — n, n 6 N. By 
construction, v is the (unique) tangency point of n + 1 Euclidean cylinders, where v belongs 
to the non-singular component of each such cylinder. Hence, such a cylinder contributes 7r 
to the cone angle at v. Recall that each such cylinder is attached (by an isometry) to the 
quotient of the bottom of another Euclidean cylinder. It is easy to check that the contribution 
to the cone angle at v from this Euclidean cylinder (with one singular boundary component) 
is equal to n x n + tt. Therefore, 

(4.12) (f>(v) = 2(n + l)?r. 

□ 
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Remark 4.13. One may also prove Theorem |0.2 by an induction on the number of boundary 
components. However, the assertion of Proposition [2728 must be used as well as an extension 
of go to <7i over the singular curve L±. Theorem 2.20 needs to be used in order to prove 
equality of the /i-length of L\ according to both the go and the g\ metric. Overall, we found 
the proof which does not use induction conceptually more gratifying. We could also stop the 



process once a planar pair of pants is encountered, thus using Theorem 0.1 directly. Finally, 
Theorem 0.1 is of course a special case of the theorem above. Still, we maintain that this 



special case deserves its own proof. 

Remark 4.14. There is a technical difficulty in our construction if some pair of adjacent 
vertices of has the same g-value (the first occurrence is in Equation ( J2. l[ ) ) ■ One may 
generalize the definitions and the index formula to allow rectangles of area zero, as one 
solution. For a discussion of this approach and others see [23J Section 5]. Experimental 
evidence shows that when the cell decomposition is complicated enough, even when the 
conductance function is identically equal to 1 and the cells are triangles, such equality rarely 
happens. 

Remark 4.15. The existence of singular curves for g results in the fact that some rectangles are 



not embedded in the target. This is evident by the proofs of Theorem 0.1 and Theorem 0.2 



Since some of the cylinders constructed have a singular boundary component, it is clear that 
some points in different rectangles that lie on this level curve will map to the same point. 
However, this occurs only in the situation above and since this fact is not of essential interest 
to us, we will not go into more details. 




FIGURE 4.16. Level curves for a 5 connected domain 
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